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摘 要 

有 限 温 度 下 自 旋 半 经 典 系统 的 随机 动力 学 行为 通常 由 随机 Landau-Lifshitz 方程 描 
B. 本 文 在 庆 之 万 随机 微分 方程 的 框架 内 ， 推 导出 有 效 朗 之 万 方程 的 一 般 形式 ， 及 其 对 
应 的 Fokker-Planck 方程 的 表达 式 . 该 有 效 明之 万 方程 能 正确 描述 正则 系 综 下 自 旋 半 经 
典 系统 的 统计 物理 性 质 ， 并 且 在 阻尼 项 和 随机 项 消失 时 能 退化 到 自 旋 半 经 典 运 动 方程 ， 
因此 是 随机 Landau-Lifshitz 方程 的 一 种 推广 . 在 笛 卡 尔 坐 标 系 和 球 坐 标 系 中 ， 分 别 给 出 
有 效 朗 之 万 方程 的 一 般 形 式 和 对 应 的 Fokker-Planck 方程 的 显 式 表达 式 .在 球 坐 标 系 中 ， 
讨论 了 上 朗 之 万 方程 中 的 纵 场 效应 ， 并 从 方程 采取 的 形式 中 给 出 是 否 包含 纵 场 效 应 的 判断 
依据 . 最 后 ， 有 效 朗 之 万 方程 在 一 个 单 自 旋 、 定 值 外 人 磁场 的 体系 中 进行 应 用 . 对 方程 采 
取 特 定 的 形式 进行 简便 的 求解 ， 并 成 功 得 到 玻 尔 兹 曼 稳 定 分 布 ， 该 结果 也 检验 了 有 效 朗 
之 万 方程 的 准确 性 . 
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自 旋 是 粒子 的 一 种 内 让 性 质 ， 是 量子 力学 导出 的 结果 . 对 磁性 物质 和 磁性 体系 ， 自 


旋 是 必 不 可 少 的 研究 要 素 . 无 论 是 理论 分 析 、 实 验 观 测 ， 还 是 在 计算 机 技术 的 高 速 发 展 


中 兴起 的 计算 模拟 ,研究 磁性 都 需要 考虑 自 旋 效 应 .本 文 并 不 从 量子 力学 的 角度 研究 磁 


性 体系 中 的 自 旋 效 应 和 磁 学 性 质 ， 而 是 以 统计 物理 的 观点 考察 自 旋 系统 . 自 旋 并 没有 经 
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典 对 应 , 尽管 它 与 经 典 力学 中 的 角 动 量 有 相似 之 处 . 因此 , 本 文 在 半 经 典 近似 的 框架 内 ， 


考虑 自 旋 半 经 典 系统 的 统计 物理 性 质 . 在 该 框架 内 ， 自 旋 以 半 经 典 力学 量 来 处 理 . 对 于 


自 旋 半 经 典 系统 , 本 文 研 究 的 统计 分 布 是 正则 系 综 下 的 经 典 玻 尔 兹 曼 分 布 . 主要 的 工具 ， 
则 是 自 旋 半 经 典 系统 的 朗 之 万 随机 微分 方程 . 


以 3 维 疝 量 S 表 示 半 经 典 系 统 中 的 自 旋 变量 .系统 的 半 经 典 运动 方程 为 3] 


这 里 我 们 采用 笛 卡 尔 坐 


-| s, | 的 标记 . 方程 中 (S) 是 系统 的 哈密 顿 量 ， 


-22 e fies 的 有 效 场 向 量 ，y 是 磁 旋 比 ，x 是 向 量 又 乘 运算 . 该 方程 有 多 种 推导 和 理 


解 的 方法 ， 如 基于 量子 力学 基本 原理 的 推导 书 3， 采 用 泊 松 括号 半 经 典 近似 的 推导 所 ， 以 


及 类 比 经 典 力学 中 角 动 量 和 力矩 之 间 的 关系 、 从 而 采用 磁力 矩 的 描述 中 等 等 ， 这 里 不 具 


体 展 开 讨 论 . 在 有 限 温度 下 , 描述 自 旋 系统 随机 运动 最 常用 的 工具 是 随机 Landau-Lifshitz 


方程 ea 


dS 031 OH 
mE (sse (- RJ us S xh(t). (2) 


个 朗 之 万 随机 微分 方程 ， 在 半 经 典 运动 方程 的 基础 上 加 入 了 阻尼 项 和 随机 项 . 其 


Hy eRe Al, h(t) 是 高 斯 白 噪声 随机 向 量 ， 满 足 
(n(D)=0, (r(e) (e) =8 8l) 3) 


(h(t) =0, (n(e)h" (0) -o(r- n1, (4) 


1 是 单位 矩阵 ， 随 机 过 程 前 面 的 系数 Ku, ， 一 般 需 要 满足 涨 落 - 耗 散 关系 pat 
MESURA 28 — exp[ —BH (S)] 是 该 随机 微分 方程 的 稳定 解 . GX B= VT s ky FE 


2 


WIRES HM, TRUE, Z 是 归 一 化 常数 . 随机 Landau-Lifshitz 方程 是 一 个 横 场 的 运 


动 方程 ， 方 程 本 身 维 持 自 旋 的 模 不 变 ， 因 此 不 包含 纵 场 效 应 ， 当 不 加 入 随机 项 ， 或 者 系 
统 处 在 零 温 下 随机 项 消失 , 剩 下 半 经 典 运 动 方程 和 阻尼 项 两 部 分 , 则 随机 Landau-Lifshitz 


方程 变 成 决定 性 的 Landau-Lifshitz-Gilbert 方程 tt 12151, 


本 文 的 主要 研究 对 象 是 用 于 描述 自 旋 羊 经 典 变量 随机 运动 的 有 效 朗 之 万 方程 . 文中 


将 具体 讨论 该 有 效 朗 之 万 方程 的 一 般 形 式 , 及 其 对 应 的 Fokker-Planck 方程 的 稳定 解 . 由 


此 得 到 的 结果 ， 可 看 作 是 对 随机 Landau-Lifshitz 方程 的 推广 和 补充 . 对 有 效 朗 之 万 方程 


一 般 形 式 的 考虑 ， 以 满足 两 点 要 求 为 基础 : 1. 以 经 典 的 玻 尔 兹 曼 分 布 为 稳定 分 布 ， 从 而 


可 以 正确 描述 自 旋 体系 在 正则 系 综 下 的 统计 性 质 ，2. 阻尼 和 随机 涨 落 这 两 项 如 果 消 失 ， 


朗 之 万 方程 将 退化 到 自 旋 半 经 典 运动 方程 (1) 式 . 本 文 的 第 2 节 , 将 在 笛 卡 尔 坐 标 系 中 根 


据 以 上 两 点 要 求 ， 推 导 有 效 开 之 万 方程 一 般 形式 的 表达 式 . 第 3 节 将 把 第 2 节 得 到 的 结 


果 ， 表 达成 球 坐 标的 形式 ， 得 到 球 坐 标 系 下 的 有 效 表 之 万 方程 和 对 应 的 Fokker-Planck 


方程 . 在 球 坐标 系 形式 下 , 可 以 比较 简便 地 讨论 和 判断 朗 之 万 方程 中 的 纵 场 效 应 . 此 外 ， 


文中 以 一 个 具体 的 简单 系统 为 例子 ， 展 示 了 选取 有 效 表 之 万 方程 的 特定 形式 并 加 以 求解 


的 过 程 . 选取 适当 的 形式 令 方 程 的 求解 、 对 系统 统计 分 布 的 分 析 变 得 简便 . 最 后 第 ACH 


是 本 文 的 总 结 . 


2 自 旋 半 经 典 朗 之 万 方程 和 Fokker-Planck 方程 
考虑 一 个 由 朗 之 万 方程 描述 的 系统 ， 以 S 表 示 系 统 的 变量 . 在 我 们 研究 的 具体 例子 


中 ，S 则 为 自 旋 向 量 .， 朗 之 万 方程 一 般 采 用 以 下 形式 
dS _ 
dt 

其 中 a(S) 是 与 $ 维 数 相 同 的 向 量 ，B(S) 是 随机 过 程 h(?) 前 的 系数 矩阵 ，h(7) 是 高 斯 随 


a(S)+B(S)h(?), (5) 


机 向 量 ， 与 (2) 式 中 的 随机 过 程 性 质 相同 ， 满 足 (3) 式 和 (4) 式 . h(t) 可 看 成 维 纳 过 程 对 时 
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间 的 形式 导数 . 朗 之 万 方程 (5) 式 , 由 决定 性 部 分 a(S) 与 随机 过 程 部 分 B(S)h(?) 组 成 . 如 


果 系 数 B(S) 由 系统 变量 $ 决 定 ， 则 该 随机 部 分 称 为 乘 性 (multiplicative) 噪声 ， 否 则 称 


为 加 性 Cadditive) 噪声 (例如 B 为 常数 矩阵 ) . 系统 遵循 朗 之 万 方程 做 随机 运动 时 ， 其 


售 时 概率 密度 分 布 p(S,7) 的 时 间 演化 规律 ， 通 常 由 著名 的 Fokker-Planck 方 程 来 描述 ， 对 


朗 之 万 方程 采取 不 同 的 随机 积分 方式 ， 相 应 地 会 导出 不 同形 式 的 Fokker-Planck 方 程 ， 最 


常见 的 有 It 随机 积分 形式 和 Stratonovich 随 机 积分 形式 . 这 两 种 积分 形式 在 处 理 


noise-induced drift 项 时 可 能 产生 差异 ， 由 此 导出 不 同 表 达 式 的 Fokker-Planck 方 程 . 对 于 


乘 性 噪声 的 随机 过 程 ， 这 个 差异 不 为 0， 对 于 加 性 噪声 这 个 差异 则 不 存在 . 由 于 通常 在 


物理 应 用 中 Stratonovich 随 机 积分 更 被 青睐 , 这 里 我 们 采取 Stratonovich 随 机 积分 所 导出 的 


形式 ， 该 形式 的 Fokker-Planck 方 程 为 0619] 


这 里 8, 指 S 的 第 ;个 分 量 ，B,(S) 指 B(S) 第 ; 行 第 大 列 的 矩阵 元 ， 容 易 看 出 ， 对 于 加 性 


噪声 的 随机 过 程 ， 方 程 中 右 端 的 第 2 项 为 0， 剩 余部 分 是 Ito 随机 积分 与 Stratonovich 随 


机 积分 的 共同 结果 . 把 方程 的 右 端 写成 时 间 演 化 算 符 Z 的 形式 ， 则 

0 E 8B,(S)| 15. TA) 
a | A sy 
概率 密度 分 布 函 数 p(S,t) 从 0 时 刻 的 初始 分 布 〈 例 如 p(S,0)=6(S-Su) ) 出 发 ， 按 照 
Fokker-Planck 方程 (0) 式 进行 演化 ， 最 终 到 达 长 时 极限 的 稳 态 p(S,w) ， 也 就 是 


Fokker-Planck 方程 的 稳定 解 . 通常 称 p(S,%m) 为 稳定 分 布 或 不 变 分 布 ， 本 文 以 pu(S) 标 


记 稳定 分 布 , 下 标 st 表示 稳定 Cstationary) . 显然 p, (S) 对 时 间 的 导数 是 0, 因此 把 o, (S) 


代入 Fokker-Planck 方程 (6) 式 的 右 端 结果 是 0， 即 Lp, (S)=0. 


本 文 研究 的 系统 是 自 旋 半 经 典 系统 . 本 文 的 研究 框架 是 探讨 适用 于 自 旋 半 经 典 系统 
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的 明之 万 方程 ， 使 其 以 经 典 玻 尔 兹 曼 分 布 为 稳定 解 ， 并 在 阻尼 项 和 随机 项 消失 时 能 回 到 


自 旋 半 经 典 运动 方程 (1) 式 . 由 此 出 发 针对 自 旋 系统 这 一 具体 例子 ， 朗 之 万 方程 的 形式 


可 以 定 为 
dS 031 


半 经 典 运动 方程 和 阻尼 项 一 起 构成 方程 的 决定 性 部 分 ，y, 是 阻尼 系数 . 随机 部 分 仍然 写 


成 系数 矩阵 B(S) 与 高 斯 白 噪声 随机 向 量 h(t) 的 乘积 ， 为 方便 起 见 我 们 加 入 了 一 个 常 关 


系数 Vs ， 类 似 随机 Landau-Lifshitz 方程 2) 式 ， 朗 之 万 方程 采取 这 种 形式 ， 显 然 满足 


在 阻尼 和 随机 两 部 分 消失 时 回 到 半 经 典 运动 方程 (1) 式 的 要 求 . 再 分 析 其 对 应 的 


Fokker-Planck 方程 . 根据 (6) 式 可 以 直接 写 出 


ð Q OH 
£o)» Es (- Sois.) |-Z ras) 0660) 
(9) 
1 8 OB, (S 1 à | LS Op(S,t 
Ju OZE a TETEN (s)s' (s) £e 
定义 向 量 
».(s) 
OB , (S 
TIR nin A ). (10) 
b, (S) r 
则 (9) 式 右 端 第 3 项 显然 有 
3B, (S oB, (S 
Ys. (52720) -y p, (s)y - yp, (s) (s)-[m(s)n(S)]. qn 
于 是 可 以 把 (9) 式 写成 更 紧凑 的 形式 
Q Q 031 
Pu )--E be (- oe J 
(12) 


系统 的 含 时 概率 密度 函数 p (S, t) 满足 这 个 时 间 演 化 方程 .对 p (St) EREKTI F H 
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稳定 解 Py (S) FE BUR 


Ii 
us 
NE 
or 
© 
yo 
a 


式 代 入 右 端 
e L l sm 
as |’ B as 
a |l 1 OH va) 
E ra (s) SB(S)b(s)+ S48(s)B" (S) Ho) -0 
为 0， 一 个 自然 的 选择 方 


其 中 的 第 1 项 ， 由 问 量 代数 分 析 不 难得 出 为 0. 第 2 ERI 


式 是 
a(S) B(8)B"() - on | J+ BB(S)b(S) (14) 
以 及 
xs 15 
Hs B (15) 


曼 分 布 为 该 Fokker-Planck 方程 的 稳定 


(15) 式 是 涨 落 - 耗 散 关 系 . 这 样 ， 则 满足 以 玻 尔 效 
解 . 这 时 上 朗 之 万 方程 的 表达 式 (9) 式 可 具体 写 为 


dS _ 031 r 031 | 1 27s 
n -x (- Jes. ns) s)|- a jns») - 5 B Ohl). (16) 
同时 相对 应 的 Fokker-Planck 方程 的 表达 式 可 以 整理 成 
6 6 031 
2»&0-- E us (- S ots t) 
2 


待定 ， 适 当地 选取 矩阵 


(16) 式 和 (17) 式 ， 是 一 个 有 效 的 统一 框架 ， 表 达 式 中 的 矩阵 B(S) 待 定 . 


均 可 以 得 到 描述 自 旋 半 经 典 系统 在 正则 系 综 下 统计 行为 的 有 效 朗 之 万 方程 


B(S), 


FEAS SCTE BER IR Ag Ps PF Bl A ER. 


UFR 26 HEM B(S) 为 


容易 验证 ， 


17) 


B(S)=| S 0 -S, |, (18) 
-S, S, 0 
则 对 任意 向 量 y 有 
B(S)y=Sxy, B” (S)y - -B(S)y - -Sxy, (19) 


同时 有 bf(S)= 0， 朗 之 万 方程 (16) 式 这 时 变 成 


B ex aslan El Prog 
TE | x y, S s | = J^ r h(t). (20) 


这 正 是 随机 Landau-Lifshitz 方程 (2) 式 . 可 见 随机 Landau-Lifshitz 方程 可 以 看 作 是 这 里 得 


到 的 有 效 朗 之 万 方程 的 一 种 具体 形式 , 本 文 把 (18) 式 中 对 B(S) 的 选取 标记 为 Bu (S) . 同 


时 可 以 把 Bu (S) 代 入 (17) 式 ， 直 接 写 出 其 Fokker-Planck 方程 的 显 式 表达 式 


(21) 
= x x _3H Vs x , 20e(S.1) 
yS s | js， S s = | 


之 前 的 很 多 研究 工作 ， 己 经 对 该 表达 式 进行 了 详细 的 推导 和 分 析 % ”31， 这 里 可 把 它 看 


作 一 个 统一 框架 的 一 种 具体 的 形式 . 


3 有 效 朗 之 万 方程 的 球 坐 标 形式 
分 析 随 机 Landau-Lifshitz 方程 (20) 式 ， 不 难看 出 ， 自 旋 变 量 的 模 |S| 在 方程 中 保持 不 


变 ， 因 为 $ 随 时 间 的 变 元 在 方向 上 与 $ 午 直 ，| 人 | 在 时 间 演 化 中 不 发 生 改变 .于 是 求解 随 
机 Landau-Lifshitz 方程 ， 将 得 到 S$ 在 一 个 固定 半径 的 球面 上 的 概率 分 布 . 这 表明 随机 


Landau-Lifshitz 方程 只 包含 横 场 作用 ， 而 没有 纵 场 效应 ， 对 于 球面 上 的 概率 分 布 ， 一 个 


自然 的 想法 是 把 系统 变换 到 球 坐 标 系 下 ,分 析 随 机 微分 方程 的 形式 及 其 对 应 的 概率 密度 


函数 . 本 贡 将 推导 球 坐 标 系 中 有 效 朗 之 万 方程 的 一 般 形 式 ， 得 到 其 对 应 的 Fokker-Planck 


方程 ， 并 判断 朗 之 万 方程 中 是 否 体现 纵 场 效 应 . 
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3.1 球 坐 标 系 中 的 有 效 朗 之 万 方程 和 Fokker-Planck 方程 


轴 的 夹 角 (0<0<zx)， 


S 
tL. PRE Syn 
9 


为 


2 是 $ 在 xz 平面 上 的 投影 向量 


S 
E 标 表示 为 | 9 |, 
9 


a 5 x 4 


RH (0€ o «2z) . 为 方便 


=| 9 |, FES sph 表示 spherical. BRA 


“ 标 变量 


& HBA 


S, - Ssin0cosg, S, - SsinOsing, S, = Scos0. 


[之 


Es SE BS) JE RT EEIE PE CORSA EE) E E 


| OS, AS, aS, 
OS 00 ô . P 
sindcosg Scos@cosg -—SsinOsingo 
OS, OS, OS, uM , 
C= — ——-|-!sinOsino Scos@sing Ssin@cosg |, 
OS 00 Od 
cos 0 —S sind 0 
os, os, OS, 
| OS 00 O@ | 
和 
OS. 0S, OS, sin cos o sin sino cos Ó 
Cre og, oe, We = : cos 0 cos o 1 kino "T i 
OS. OS, os, S S S 
0p Op 09 一 sing —— cos 9 0 
os, oS, os. | SsinO Ssin@ 4l 


IC| = S^ sind 4 


其 中 


利用 变换 矩阵 ， 可 以 方便 


Pa (Son) = 


Z 
定义 有 效 的 球 坐 标 哈密 


1 


E 标 变换 的 雅 可 比 行列 式 . TERA 


8 


Abn Soon » KRESS 


p| -BH (s) |? sind =—exp| - -PH (S sa) 


)=H(S)—n(s" sin). 


分 布 的 形式 为 


地 把 上 一 节 得 到 的 朗 之 万 方程 (16) 式 转换 到 球 坐 标 系 


其 中 5S 是 S$ 的 模 ，09 是 $ 与 笛 卡 尔 坐 标 系 中 z 


4 标 变 量 间 的 关系 


(22) 


(23) 


(24) 


(25) 


(26) 


Boy oi Lucus. ( aH rczf 020) 1 EA 
"xx 274€ s | E Jes oet T Jets). REED 


这 里 定义 D=C-B. 逐 项 分 析 上 式 的 右 端 . 首先 对 应 随机 Landau-Lifshitz 方程 中 的 Bi 


[ 见 (18) 式 ]， 有 


0 0 0 
D,-CB,-| sing —cosp O0. (28) 

cotücosp cot@sing -l1 
容易 验证 


0 0 0 
1 os eS, Os 
D,, - — 一 二 一 一 一 一 一 二 |. (29) 
Ssin@| 60 Op Op 


于 是 (27) 式 右 端 第 1 项 为 


| 
本 加 | 


(30) 
1 OS, oS, 60S, | A 1 | 031 


ano oa op Sell S|)’ Sna om 


对 第 2 项 ， 首 先 有 


as, oS, as 
ôS ôS as 
‘(-)- ôS, OS, OS, HE | 031 GD 
os 00 00 00 əs ôS yn 
ôS, ôS, 68, 
Lep 6p 0g | 


再 分 析 向 量 b(S) .在 (10) 式 中 定义 了 b(S) 的 元 素 ， 这 里 用 球 4 


6 标 变量 来 表示 b(S) 


OB, OS, 
b, (S) QU eI EJ -FT Ec) 


sph,i 


dzee) aD 
LOG = LE a LG TALES AG 
S ahi Sani J S ahi (32) 
aD, 0C, Sou, aD, oC , 
22215 CTUM 2D; 25» os S, “Lig +È Di F os, 
2/S 
= 5 Pu yr cot Ó |, 
最 后 一 步 用 D; 表示 D7 AN kIT, HAA M23) sR E BI ER 
C, 
25s 
2/8 
CCB tO (33) 
- es. —| co i 
j 0 
eC, 
748. 
由 (32) 式 显然 可 以 把 向 量 b(S) 表达 成 更 紧凑 的 形式 
2/5 
b(S)=d(S,,,)+D"| cotà |, (34) 
0 
其 中 定义 向 量 d(S。, ) 为 
d, (S...) 
OD, 
d(S,,,)=| d, (Saa) |. d, Ga-x os (35) 
a, (S) = 


可 以 看 出 定义 方式 与 b(S) 在 笛 卡 尔 坐 标 系 中 的 形式 相似 .现在 把 (30)、(31) 和 (34) 式 代 


入 (27)， 得 到 


2/8 
dS 
wh zd +y; 3 DD” E 7t ++pd(S,,,)}+ E ph (1) 
dt Ssin@| dg OS.) P A B p B 


(36) 


031. 
-y E «y, | DD' -Sm |++a(s,, )|+ E p(t). 
SsinO| 09 OS... ) P n p 


这 是 有 效 朗 之 万 方程 (16) 式 的 球 坐 标 形 式 . AREE BH, 现在 转化 为 对 了 D 的 选取 . 对 


应 随机 Landau-Lifshitz 方程 的 取 法 ， 是 (28) 式 中 的 D,，. 这 时 含 时 概率 密度 函数 


Pon (Sonst) 的 时 间 演化 ， 则 满足 球 坐标 系 的 Fokker-Planck 方程 ， 与 第 2 节 的 分 析 完 全 


类 似 地 ， 可 以 得 到 (36) 式 对 应 的 球 坐 标 系 Fokker-Planck 方程 


" S f 
as... P Sind dg Pul sph ) 
OH (37) 


3H Op (Seon ot 
+7sDD* | -——= Pon (Sat) 22 DD" Pal : ) i 
Son mo Pp OS 


不 难 检验 ， 球 坐标 的 玻 尔 兹 曼 分 布 (25) 式 是 该 Fokker-Planck 方程 的 稳定 解 . 把 (25) 式 代 
入 (37) 式 右 端 ， 得 到 


0 
ð 1 8e 031, 0H, 
as, | BZ 6g nd - as “Jp wu (S.,.) - 7, DD' - ag. |^ on (Son) G9 
sph sph sph 
1 de 
| \ pz 00 ] 


结果 显然 是 0. 这 表明 该 朗 之 万 方程 在 球 坐 标 系 中 的 有 效 性 . (36) 式 和 (37) 式 ， 是 有 效 朗 
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之 万 方程 及 其 Fokker-Planck 方程 的 一 般 形式 (16) 和 (17)， 在 球 坐 标 系 中 的 表现 形式 .这 


是 本 文 在 球 坐 标 系 中 得 到 的 主要 结果 . 


3.2 纵 场 效应 的 讨论 


从 球 坐标 形式 的 有 效 朗 之 万 方程 (36) 式 可 以 直接 看 出 ， 如 果 只 有 横 场 作用 ， 保 持 模 


5S 守恒 , 则 必须 满足 矩阵 D 的 第 1 行 是 零 向 量 .随机 Landau-Lifshitz 方 程 的 Dii [ 见 (28) 式 ]， 


是 满足 该 条 件 的 一 个 例子 . 把 随机 Landau-Lifshitz 方程 对 应 的 球 坐 标 形式 显 式 地 写 出 


0 0 
0 0 0 
dS 3H 
may l 091 +y; ERG T + A sing -cosp 0 |h(t).(39) 
dt Ssin@| OQ 00 p 10 ‘ei 1 
cotÜcoso cot@sing 一 
031 1 0H, ? á 
00 sin? 0 0Q 


显然 二 =0， S 是 常数 ， 方 程 转化 成 两 个 角度 变量 的 随机 微分 方程 


T 5 a «| Ti "im [^5 [sin oh, (1) -coseh; (£)] 


(40) 
d 1 0H 1 0H, [2 
=a aa Ara ary mee ey 22 E | cot 8 cos gh, (1) + cot sin gh, (t) — h (1) . 


相应 地 ， 概 率 密度 函数 p,,, (9,9,1) 只 包含 两 个 角度 变量 和 时 间 变 量 ，Fokker-Planck 方程 


可 整理 得 到 
0 0 0 
== E "X : 031 Pan tY - $31, p E Peon 
0t " 0S. | Ssind| dg | " ^ 00 ™ g 00 
u OH 1 OH on 1 Pon 
| 00 sin^9 dp sin’ 0g 
o 1 OH 031, K S 0p, sph 
= + gem 41 
2l "Ssmo Oo 9" /* sp "-" 5 oo oo 


4 Q 1 OH 1 OFT mE + Vs 1 OP son 
spl” Sand od a6 to ^" Bane bo | 


这 个 结果 已 经 被 很 多 工作 分 析 和 讨论 过 P2 264), BEDS 1 行 如 果 不 是 零 向 量 ， 则 上 朗 
之 万 方程 中 包含 纵 场 效应 ， 模 S 不 能 保持 不 变 .在 文献 [35] 中 ， 作 者 提出 了 一 个 包含 纵 
场 涨 落 的 朗 之 万 方程 ， 并 说 明 随 机 Landau-Lifshitz 方程 是 该 方程 在 以 S 为 半径 的 球面 上 


的 投影 . 在 本 文 的 框架 中 ， 该 方程 是 (160) 式 对 矩阵 B 简单 选取 为 B=1 得 到 的 结果 ,或 者 


等 价 地 在 球 坐 标 系 中 选取 D=C  . 这 样 选取 的 矩阵 D 的 第 1 行 显然 不 是 零 向 量 ， 因 此 
方程 中 包含 纵 场 效应 . 除了 朗 之 万 方程 以 外 ， 纵 场 效 应 的 研究 也 可 见于 蒙特 卡 洛 模拟 的 
HRT PLES, 


如 果 只 考虑 横 场 作用 ， 有 效 朗 之 万 方程 也 可 以 采取 随机 Landau-Lifshitz 方程 以 外 的 


形式 . 在 文献 [32] 中 ， 作 者 采用 了 


00 0 
D=/0 1 0 (42) 
0 0 E 
L sin Ó | 
而 得 到 
031, 
dO | l OH — : sph + 27s, (1) 
dt SsinÓ OQ 00 Pp (43) 
do 1 OH 1 031 on 275 1 
—=-y - he 7 二 | 一 一 === h, (1). 
dt Ssin@ 00 sin 0 OQ p sin 


0 


可 以 看 出 ， 该 形式 的 朗 之 万 方程 与 随机 Landau-Lifshitz FEOIR LL, PAP fü EAE 


jw 


在 动力 学 演化 时 都 发 生 耦 合 ， 即 运动 方程 不 独立 . 主要 区 别 则 在 随机 部 分 ， 该 方程 只 需 


要 一 个 2 维 的 高 斯 随机 向 量 场 ， 随 机 Landau-Lifshitz 方程 则 需要 完整 的 3 维 高 斯 随机 过 


Ke. 对 该 方程 具体 的 分 析 和 数值 计算 ， 这 里 不 详细 讨论 . 


3.3 简单 的 应 用 
本 文 已 经 得 到 适用 于 自 旋 半 经 典 系统 的 有 效 朗 之 万 方程 的 一 般 形式 , 在 笛 卡尔 坐标 


系 和 球 坐 标 系 中 分 别 展示 了 有 效 朗 之 万 方程 及 其 对 应 的 Fokker-Planck 方程 的 表达 式 .对 


于 具体 的 系统 、 有 具体 的 哈密 顿 量 ， 选 取 特 定 的 形式 可 能 会 令 方程 的 求解 变 得 简便 .本 贡 


将 讨论 一 种 较为 简单 的 体系 ， 对 该 体系 选取 有 效 朗 之 万 方程 的 一 种 形式 并 加 以 求解 . 


考虑 一 个 单 自 旋 的 系统 ， 处 于 定 值 外 磁场 了 再 当中 ,系统 的 哈密 顿 量 为 自 旋 与 外 场 的 


相互 作用 


H (S) - -S-H. (44) 


由 于 于 是 定 值 向 量 ， 为 方便 起 见 可 把 再 的 方向 设 为 z 轴 正 方向 ， 则 哈密 顿 量 为 


H -—SH cos0, AUR p. 球 坐 标 哈密 顿 量 有 简单 的 形式 


H. (Sp) = 5H cos0 -—In( 5" sin 0). (45) 


Sp! 


代入 有 效 朗 之 万 方程 的 一 般 形式 (36) 式 ， 得 到 


Le 


BS 

dS ° 1 1 2 

=) 0 [+y,] DD"| -SHsind+—cot@ |-—Dd(S,,) -,|—5Dh(/). (46) 
dt B B B 


为 确保 系统 有 合法 的 玻 尔 效 曼 分 布 ， 我 们 不 引入 纵 场 效应 ， 于 是 8 是 常数 ， 同 时 系统 的 


BUR BK 89 d PASE. OEE DD ， 我 们 选取 非常 简单 的 形式 


0 0 0 
0 1 0l (47) 
0 0 1 


由 此 得 到 的 朗 之 万 方程 的 具体 表达 式 为 


do 7 2 
um -7s [SH sina on 中 FA (t) 


(48) 
8. opua [s 
"x vH + B h, (7). 
可 以 看 到 ， 两 个 角度 变量 的 运动 方程 是 独立 的 朗 之 万 方程 . 对 这 两 个 方程 分 别 求 解 ， 就 


得 到 g 和 gp 各 自 的 概率 分 布 . 


首先 分 析 6 的 运动 方程 


031. 2 
de =-y,| SH NE 4 DA (t) ——ys Ur SAY (t). (49) 
dt p B 00 


一 个 过 阻尼 形式 (overdamped) 的 天 之 万 方程 . 显然 ， 这 个 方程 有 加 性 噪声 的 随机 


-— 
过 程 ， 随 机 项 前 面 的 系数 和 阻尼 因子 满足 涨 落 - 耗 散 关系 .利用 第 2 节 的 分 析 方法 ， 容 
易 知道 该 朗 之 万 方程 得 到 的 稳定 分 布 是 g 的 玻 尔 兹 曼 分 布 却 exp[-B3fw (0) ] .再 看 p 的 


ap __ EA 
"ug yH + 3 h, (1). (50) 


这 个 方程 的 形式 解 可 以 直接 写 出 


运动 方程 


9(1)=0,— yes IP Qs (51) 


Feb gy EMI EE. NIRS, o(r) BROKER, BEY, -yHt THN 


ons yar] e nsns nas) =f fn Ca 


(52) 
bei 
BU 
[B | [o7 (o rH)] 
p(o.t) - "pm «| ag (53) 


如 果 g 是 实数 域 上 的 无 界 变 量 ， 则 上 式 的 售 时 概率 密度 在 长 时 极限 下 只 有 平庸 解 0， 不 


存在 稳定 分 布 . 但 在 球 坐 标 系 中 , gp 被 限定 在 [0,2x) 区 间 上 , 朗 之 万 方程 的 形式 解 (51) 式 


中 应 包含 周期 边界 条 件 . 显然 ， 在 周期 边界 条 件 的 作用 下 ， 随 机 变量 g(1) 的 均值 也 在 


a a 


间 上 的 等 概率 密度 ， 从 而 ， p 的 稳定 分 布 是 [0,27) 区 间 上 的 均匀 分 布 二 . 于是， 系统 


在 球 坐标 系 中 总 的 稳定 分 布 可 以 得 到 


Pon (8,9) = oe | -LH (8) |= ze [ BSH cos 0|? sin 6. (54) 


这 正 是 球 坐 标 系 中 的 玻 尔 兹 曼 分 布 . 


本 节 对 具体 的 体系 进行 分 析 ， 选取 有 效 朗 之 万 方程 的 一 个 特定 的 形式 ， 从 而 可 以 相 


对 简便 地 求解 ， 以 得 到 稳定 的 玻 尔 兹 曼 分布 . 这 是 对 文中 得 到 的 有 效 朗 之 万 方程 一 般 形 


式 的 一 个 简单 的 应 用 . 此 外 ， 横 场 自 旋 半 经 典 朗 之 万 方程 在 随机 Landau-Lifshitz 方程 之 


外 有 更 多 可 能 的 形式 ， 本 贡 提 供 了 一 个 具体 的 例子 . 


4 结论 


本 文 以 正则 系 综 下 的 经 典 玻 尔 效 曼 分 布 为 目标 ， 从 明之 万 随机 微分 方程 出 发 ， 推 导 


出 自 旋 半 经 典 系 统 有 效 朗 之 万 方程 的 一 般 形式 ， 以 及 对 应 的 Fokker-Planck 方程 .该 有 


效 姑 之 万 方程 以 玻 尔 兹 曼 分布 为 稳定 分 布 , 因此 可 以 正确 描述 自 旋 体系 在 正则 系 综 下 的 


统计 物理 性 质 . 同时 方程 中 阻尼 项 和 随机 项 消失 时 能 退回 到 自 旋 半 经 典 运动 方程 ， 因 此 


也 能 包含 自 旋 系统 的 半 经 典 运 动 模式 . 这 个 结果 ， 是 对 随机 Landau-Lifshitz 方程 的 推广 


和 补充 . 在 球 坐 标 系 中 , 方程 的 形式 可 以 简便 地 判断 是 否 包含 纵 场 效应 . 在 一 个 单 自 旋 、 
定 值 外 磁场 的 具体 体系 中 ， 成 功 地 应 用 有 效 朗 之 万 方程 进行 统计 分 布 的 分 机， 检验 了 方 


程 的 准确 性 ， 同 时 也 为 横 场 自 旋 朗 之 万 方程 在 随机 Landau-Lifshitz 方程 之 外 更 多 的 形式 


上 的 选择 ， 提 供 了 具体 的 例子 . 本 文 的 工作 ， 可 以 拓宽 朗 之 万 随机 微分 方程 作为 一 种 理 


论 和 计算 工具 的 应 用 范围 , 同时 也 为 自 旋 半 经 典 系统 的 动力 学 和 统计 力学 研究 提供 理论 


工具 层面 的 参考 . 
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Abstract 

The stochastic dynamics of spin semiclassical system at finite temperature is usually 
described by stochastic Landau-Lifshitz equation. In this work, the stochastic differential 
equation for spin semiclassical system is studied. The generalized formulation of effective 
Langevin equation and the corresponding Fokker-Planck equation are derived. The obtained 
effective Langevin equation offers an accurate description of the distribution in the canonical 
ensemble for spin semiclassical system. When the damping term and the stochastic term 
vanish, the effective Langevin equation reduces to the semiclassical equation of motion for 
spin system. Hence, the effective Langevin equation can be seen as a generalization of the 
stochastic Landau-Lifshitz equation. The explicit expressions for the effective Langevin 
equation and the corresponding Fokker-Planck equation are shown in both Cartesian and 
Spherical coordinates. It is demonstrated that, the longitudinal effect can be easily illustrated 
from the expressions in Spherical coordinates. The effective Langevin equation is applied to 
the simple system of a single spin in a constant magnetic field. In choosing an appropriate 


form, the Langevin equation can be easily solved and the stationary Boltzmann distribution 
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can be obtained. The correctness of the Langevin approach to the spin semiclassical system is 


thus confirmed. 
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